Signature of smooth transition from diabatic to adiabatic states 
in heavy-ion fusion reactions at deep subbarrier energies 
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We propose a novel extension of the standard coupled-channels framework for heavy-ion reactions in order 
to analyze fusion reactions at deep subbarrier incident energies. This extension simulates a smooth transition 
between the diabatic two-body and the adiabatic one-body states. To this end, we damp gradually the off- 
diagonal part of the coupling potential, for which the position of the onset of the damping varies for each 
eigen channel. We show that this model accounts well for the steep falloff of the fusion cross sections for the 
i6o+208pt,^ «Ni+^''Ni, and '^'Ni+^^^Ni reactions. 

PACS numbers: 25.60.Pj, 24.10.Eq, 25.70.Jj,25.70.-z 



Heavy-ion fusion reactions at low incident energies provide 
a good opportunity to study fundamental features of the tun- 
neling phenomena in many-particle systems. A potential bar- 
rier, called the Coulomb barrier, is formed because of a strong 
cancellation between the repulsive Coulomb interaction and 
an attractive nuclear interaction. In particular, the potential 
tunneling at incident energies below the Coulomb barrier is 
refeiTed to as the subbarrier fusion reaction. One prominent 
feature of the subbaiTier fusion reactions is the large enhance- 
ment of fusion cross sections, as compared to a prediction of 
the simple potential tunneling 1 1]. This enhancement has been 
attributed to the coupling of the relative motion between the 
colliding nuclei to several intrinsic degrees of freedom, such 
as a collective vibration of the target and/or projectile nuclei. 
The coupled-channels (CC) approach, based on this picture, 
has been successful in accounting for the subbarrier enhance- 
ment 1 2]. 

Because of a recent progress in experimental techniques, 
it has been possible to measure fusion cross sections down 
to deep-subbarrier incident energies |3i 0, [H |^. These data 
show a substantial reduction of fusion cross sections at deep- 
subbarrier energies from the prediction of the CC calculations 
that reproduce the experimental data at energies around the 
Coulomb baiTier, and have brought about a renewed interest 
in this field. This phenomenon, often refeiTed to as the fusion 
hindrance, shows a threshold behavior, where the data deviate 
largely from the standard CC calculations at incident energies 
below a certain threshold energy, 

A key element to understand the fusion hindrance is the 
density overlap of the colliding nuclei in the potential tunnel- 
ing process. When the incident energy is below the potential 
energy at the touching point of the colliding nuclei, Vxouch, 
the inner turning point of the potential is located inside the 
touching point, and the projectile is still in the classically for- 
bidden region when the two nuclei touch with each other (see 
Fig. 1 in Ref.[71]). In this situation, the colliding nuclei have to 
penetrate through a residual barrier with an overlapping con- 
figuration before fusion takes place. In our previous work 1*^, 
we evaluated VtoucH using several kinds of ion-ion potential, 
and systematically compared those with experimentally de- 



termined threshold energy for many systems. The obtained 
systematics shows a strong coiTelation between Vxouch and E^, 
indicating strongly that the nuclear interaction in the overlap- 
ping region plays a decisive role in the deep-subbarrier hin- 
drance. 

Three different mechanisms have been proposed so far in 
order to account for the fusion hindrance. Based on the sud- 
den picture, Mi§cu and Esbensen have investigated the effect 
of the nuclear interaction in the overlap region, in terms of a 
repulsive core due to the Pauli exclusion principle lIsllgl fToll . 
Assuming the frozen-density in the overlapping region, they 
obtained a much shallower potential pocket than the stan- 
dard one, which hinders the fusion probability for high partial 
waves. Their shallow potential reproduces well the fusion hin- 
drance. In contrast, we have proposed the adiabatic approach 
by assuming neck formations between the colliding nuclei in 
the overlap region |11]. In our model, the fusion hindrance 
originates from the tunneling of much thicker potential barrier 
characterized by the adiabatic one-body potential. This model 
achieved comparable good reproduction of the experimental 
data to the sudden model. The third mechanism, suggested 
recently by Dasgupta et al., is the quantum decoherence of 
channel wave functions caused by the coupling to the thermal 
bath ilt]. A model calculation based on this picture shows a 
possibility of the gradual occurrence of hindrance in subbar- 
rier fusion reactions [ 12J . 

Among those three mechanisms, the origin for the deep- 
subbaiTier hindrance is considerably different from each other. 
The recent precise data for the "'0-H^°^Pb fusion |5] may pro- 
vide an adequate system to discriminate among the various 
models, because the behavior of its astrophysical S-factor is 
difficult to reproduce within a simple model calculation. In the 
model of Esbensen and Mi§cu, not only the collective inelastic 
channels but also the particle transfer channel with modified 
coupling strengths are necessary for a fit to the experimental 
data [8]. In the estimation of Dasgupta et al., it was impossi- 
ble to obtain an overall fit to the experimental data from the 
above-barrier to deep-subbarrier regions with a single param- 
eter set for the nuclear potential |5]. On the other hand, the 
performance of the adiabatic model has not yet been studied 
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FIG. 1: (Color online) Fusion cross sections for the *"*Ni+*'*Ni and 
i6Q^_208pj^ systems. The solid and the dashed lines are the calculated 
result with and without the damping factor, respectively. The dot- 
dashed line is the result of no coupling with the YPE potential. 
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FIG. 2: (Color online) Astrophysical S-factor for the *'*Ni+*'*Ni and 
i6Q^208pj^ systems as a function of the incident energies. The mean- 
ing of each line is the same as in Fig. 1 except for the dot-dashed line, 
which shows the result of the sudden model taken from Refs. Llftl 
andd]. 



for this system, although the concept of the adiabatic potential 
was proved helpful in the analysis of the potential inversion 
method in the deep-subbarrier fusion fTst]. 

In this paper, we attempt to study the deep-subbarrier fusion 
for the '*0-i-^'"^Pb reaction based on the adiabatic model. Our 
previous model has a defect that the full quantum treatment for 
the two-body part suddenly switches to the semi-classical one 
for the adiabatic one-body part, which introduces arbitrariness 
for the choice of the Hamiltonian. To avoid the shortcoming, 
we develop below a full quantum mechanics where the CC 
approach in the two-body system are smoothly jointed to the 
adiabatic potential tunneling for the one-body system, result- 
ing in an overall good agreement for the ^^O+^^^Ph reaction, 
as well as for the ^^^Ni+^'^Ni and '''*Ni+^'*Ni systems. 

We employ the incoming wave boundary condition in order 
to simulate a compound nucleus formation. In order to con- 
struct an adiabatic potential model with it, we postulate the 
followings: (1) Before the target and projectile nuclei touch 
with each other, the standard CC model in the two-body sys- 
tem works well. (2) After the target and projectile overlap 
appreciably with each other, the fusion process is governed 
by a single adiabatic one-body potential where the excitation 
on the adiabatic base is neglected. (3) The transition from the 
two-body treatment to the one-body one takes place at near 
the touching configuration, where all physical quantities are 
smoothly joined. 

To this end, we adopt Yukawa-plus-exponential (YPE) po- 
tential fTT] as a basic ion-ion potential V^J", because the di- 
agonal part of this potential satisfies the conditions (l)-(3) by 
choosing a suitable neck-formed shape for the one-body sys- 
tem, as has been shown in our previous work fTT] . In addition, 
the saturation property of the nuclear matter is phenomeno- 
logically taken into account in the YPE model. It has also 
been shown that the YPE model is consistent with the po- 
tential obtained with the energy density formalism with the 
Skyrme SkM* interaction ITslflQl . 

The nuclear coupling form factor which describes excita- 



tions to the vibrational states in the two-body channel is taken 
as the derivative of potential V^' |17]. The coupling matrix 
elements are evaluated with the eigen-channel representation 
as in Eq. (24) in Ref.|18]. In order to satisfy the conditions 
(l)-(3), we employ the following form for the nuclear poten- 
tial for the eigen-channel a. 
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where Aa is the eigen value of the excitation operator. The 
most important modification from the standard CC treatment 
is the introduction of the damping factor <1>. This damping fac- 
tor represents the physical process for the gradual transition 
to the adiabatic approximation, by diminishing the strength of 
excitations to the target and/or projectile vibrational states af- 
ter the two nuclei overlap with each other. We thus choose the 
damping factor given by 

(1 r > Rd + Aa (Two-body region), 

(2) 
e r < Rd + Aa (Overlap region), 

where Rd is the spherical touching distance between the tar- 



get and projectile defined by Rd - rd(A}^^^ + Ap'), rd is the 
damping radius parameter, and Ud is the damping diffuseness 
parameter. Notice that the touching point in the damping fac- 
tor depends on Aa, that is, the strength of the excitations starts 
to decrease at the different distance in each eigen channel. 

It is slightly complicated to take into account the effect of 
the damping factor on the Coulomb coupling. When different 
multipole components are present simultaneously, the eigen 
channels, which are introduced to evaluate the nuclear cou- 
pling matrix elements, do not diagonalize the Coulomb cou- 
pling matrix. We have therefore introduced the channel in- 
dependent damping factor for the Coulomb coupling, but the 
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FIG. 3: (Color online) Logarithmic derivatives of fusion cross sec- 
tions, L(£) = dln(£c.m.o-fus)/rf£'c.m.,forthe**Ni+«''Niand '*0+2"'*Pb 
systems as a function of the incident energies. The meaning of each 
line is the same as in Fig. 1. 



effect on the fusion cross sections appeared small. For sim- 
plicitly, we therefore consider the damping factor only for the 
nuclear coupling in the calculations presented below. 

We apply our present model to the fusion reactions for the 
^'^Ni+^'^Ni and '^O+^^^Pb systems. To this end, we incorpo- 
rate the damping factor and the YPE potential in the computer 
code ccfull [18]. For the '''^Ni+^'^Ni system, the coupHng 
scheme included in the calculation, as well as the deformation 
parameters, are the same as in Ref. ||4<1 . For the '^OH-^"^Pb 
system, those are the same as in Ref. 11911 . but we include the 
coupling to the low-lying 3" phonon states and the double- 
octupole phonon excitations for both the '^O and ^"^'Pb nuclei. 
For the damping factor, we use - 1 .298 fm and - \ .05 
fm for the '"*Ni+^''Ni system, and = 1 .280 fm and = 1 .28 
fm for the "^0+^'"^Pb system. 

For the YPE model, the parameters are taken as ao = 0.68 
fm, a, = 21.33 MeV, and k, = 2.378 from FRLDM2002 ||20|l . 
In order to fit the experimental fusion cross sections, the ra- 
dius parameter ro is adjusted to be 1.205 fm and 1.202 fm 
for the ^'*Ni+'''*Ni and '^OH-^'"^Pb systems, respectively. For 
the mass asymmetric '^O-H^'^^Pb system, it is difficult to joint 
smoothly the potential energies between the two-body and the 
adiabatic one-body systems at the touching point, because the 
proton-to-neutron ratio for the one-body system differs from 
that for the target and projectile in the two-body system. To 
avoid this difficulty, we smoothly connect the potential en- 
ergy around the touching point to the liquid-drop energy of the 
compound nucleus, using the third-order polynomial function 
(see. the dashed line in Fig. 4). We do this by identifying the 
internucleus distance r with the centers-of-masses distance of 
two half spheres. The obtained potential is similar to the re- 
sult of the density-constrained time-dependent Hartree-Fock 
method ll2lll . We have checked this prescription for the mass 
symmetric ^^Ni-i-''"^Ni system, by comparing to the potential 
energy used in our previous work ifTTll . The deviation due to 
this prescription is negligibly small. 

Figure 1 shows the fusion cross sections thus obtained. The 
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FIG. 4: The adiabatic potential for the '^0-l--''*Pb system as a func- 
tion of the center-of-mass distance. The solid line is the adiabatic 
potential obtained with the damping factor. The dashed line is the 
result obtained with the conventional CC approach. The dotted line 
and the solid circle are the pontential and the touching point for the 
uncoupled case, respectively. The gray region denotes the adiabatic 
potential obtained with the potential inversion method, taken from 
Ref. 



fusion cross sections obtained with the damping factor are in 
good agreement with the experimental data for both the sys- 
tems (see the solid line). The dashed line is the YPE potentials 
without the damping factor. The dot-dashed line is the results 
of no coupling with the YPE potential. For both the systems, 
we see that drastic improvement has been achieved by taking 
into account the damping of the CC form factors. 

We also compare the astrophysical S factor representation 
of the experimental data with the calculated results, as shown 
in Fig. 2. In the calculation, the Sommerfeld parameter rf is 
shifted by 75.23 and 49.0 for the ''^Ni+'^^Ni and "'O+^'^^Pb 
systems, respectively. The S factor obtained with the damp- 
ing factor are consistent with the experimental data for both 
the systems (see the solid lines), and reproduce well the peak 
structure. 

For astrophysical interests, it is important to evaluate fusion 
cross sections at extremely low incident energies, which is dif- 
ficult to measure directly. Notice that the S factor predicted by 
our model differs considerably from that of the sudden model 
by Mi§cu and Esbensen flO], denoted by the dot-dashed line. 
For both the systems, as the incident energy decreases, their 
S factor falls off steeply below the peak of the S factor, while 
our S factor has a much weaker energy dependence. 

Figure 3 compares the logarithmic derivatives 
c/ln(£'c.m.o"fus)/<^£^c.m. of the experimental fusion cross 
section with the calculated results. It is again remarkable that 
only the result with the damping factor achieves nice repro- 
duction of the experimental data. For the ^^NiH-^'*Ni system, 
the result with the damping factor becomes saturated below 
£^c.m.=86 MeV. This behavior is similar to the experimental 
data for the '^0-i-^'"^Pb system. The measurement at further 
lower incident energies for this system will thus provide a 
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strong test for the present adiabatic model. 

Figure 4 shows the adiabatic potential of the '*0+^°^Pb sys- 
tem, that is, the lowest eigenvalue obtained by diagonalizing 
the coupling matrix at each center-of-mass distance r. The 
adiabatic potential calculated with and without the damping 
factor are denoted by the soUd and the dashed lines, respec- 
tively. The uncoupled YPE potential is also shown by the dot- 
ted line. The solid circle denotes the touching point in the 
absence of channel coupling. We see that the result obtained 
with the damping factor is much thicker than that of the con- 
ventional CC model. In this respect, it is interesting that the 
result with the damping factor is similar to that obtained with 
the potential inversion method [13], denoted by the gray re- 
gion, justifying our treatment for the damping of the CC form 
factor. 

For the average angular momentum of the compound nu- 
cleus, the results with the damping factor become saturated at 
incident energies below the threshold energy with decreasing 
incident energy, as shown in our previous works 1 1 1 , 22] . This 
result largely differs from that obtained with the sudden model 
by Mi§cu and Esbensen. Their result is strongly suppressed at 
energies below the threshold energy. It is thus interesting to 
measure the average angular momentum at deep subbarrier 
energies, in order to discriminate the two approaches. 

We have also applied our model to the ^^Ni-n^^Ni reactions 
and the results obtained are in good agreement with the exper- 
imental data. For the damping factor, we used r j = 1.3 fm and 
ad - 13 fm in order to fit the experimental data. Notice that 



the obtained damping radius parameters for the three systems 
which we study are almost the same. We emphasize that our 
model achieves an overall fit not only to the fusion cross sec- 
tions but also to the S factors and the logarithmic derivatives 
simultaneously. 

In summary, we have proposed a novel coupled-channels 
approach for heavy-ion fusion reactions by introducing the 
damping of the CC form factor inside the touch point in or- 
der to simulate the transition from the diabatic to adiabatic 
states. The important point in our present model is that the 
transition takes place at different places for each eigen chan- 
nel. By applying this model to the i^O-i-^"¥b, the ^'*Ni+^'*Ni, 
and the ^^Ni-i-^^Ni systems, we conclude that the smooth tran- 
sition from the diabatic two-body to the adiabatic one-body 
potential is responsible for the steep falloff of the fusion cross 
section. It is an interesting future study to apply the present 
model systematically to other systems and clarify the dynam- 
ics of deep subbarrier fusion reactions, and thus many-particle 
tunneling phenomena. 
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